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Abstract 

The long time asymptotics of multi-time correlation functions of relaxing quantum 
mechanical systems can be conveniently studied by means of free-products of suitable C*- 
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paper, we study the distribution law of fluctuations of temporal averages in a class of 
quantized toral automorphisms. 
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1 Introduction 



Recently, it has been suggested [BF, ABDF] that the statistics of long time asymptotics of 
multi-time correlation functions of quantum dynamical systems might be used to associate 
particular probability distribution laws with the variety of phenomena commonly going 
under the name of quantum chaos [CC]. The proposed setting is that of a free-product of 
copies of the algebra of obervables of the system, equipped with a state obtained by a 
suitable multiple time-average. 

Given an invariant state and expectations of observables at different times, with a 
same time possibly appearing more than once, one may try to express such expectations 
in terms of time-ordered correlation functions. The intrinsic complexity of the dynamics 
might render this task very hard indeed and is likely to generate a more or less unwieldy 
proliferation of commutators. We will in particular address the situation where differences 
between different times become large. In the particular case of strong time-asymptotic 
commutativity one can easily group observables at equal, largely separated, times. 

In general, the idea is to leave aside any attempt at simplifying multi-time correlation 
functions by expressing them in terms of ordered correlation functions, but rather to regard 
them as elements of a free-product algebra whose expectations are obtained by averaging 
with respect to all different times. It is plausible that the chaotic features of the dynamics, 
if any, might leave their imprinting on the statistical properties of the asymptotic state. 

Ultimately, the testing ground of the scheme of above should consist of those quantum 
systems, typically finite dimensional, that tend to classically chaotic systems when the di- 
mension of the associated Hilbert spaces increases to infinity [BV, BBTV, CC, D, HKS]. 
This involves also an appropriate rescaling of the time. In this paper, we analyse a family 
of toy models, namely the hyperbolic toral automorphisms quantized as in [BNS]. These 
models are infinite dimensional quantum systems arising via a non-commutative deforma- 
tion of the algebra of continuous functions on the two-dimensional torus. The quantum 
algebra is equipped with the tracial state and endowed with a dynamics such that the GNS 
time-evolutor coincides with the classical Koopman-von Neumann unitary operator. The 
absolute continuity of the spectrum of the quantum evolution guarantees those clustering 
properties which are ruled out in the commonly studied classically chaotic quantum systems 
which usually have discrete spectra and quasi-periodic time-behaviour. 

The mixing properties depend on the value of a certain deformation parameter 6 and 
one can distinguish different randomness conditions that lead, via a central limit theorem, 
to a variety of statistics of fiuctuations ranging from the Gaussian distribution to Wigner's 
semi-circle law. 



1 



2 Statistics of multi-time correlation functions 



We briefly resume the approach of [ABDF]. (21, 0, 0) wiU denote a discrete time dynamical 
system where 

• 21 is a unital C*-algebra 

is a discrete dynamical group of automorphisms of 21 such that 
X{t) :— ©t(X) represents the operator X evolved up to time t 

• is a ©-invariant state on 21: </> o = (/>. 

We shall consider multi-time correlation functions of the form 

t ^ 0(x(^)(t.(i))X(2)(t,(2)) . . .X(")(t.(,))) , (2.1) 

where the X^^\ti,(^j-)) are operators at times tj^(j) in 21, t = {ti,t2,---} £ and u 

maps {1,2, into Nq. Two consecutive time-indices will always be considered dif- 

ferent, otherwise, if, say — v{j + 1) — p, then we write X^^\t,y(^j^)X^^~^^\ti^(^j^i^) = 

(^X(j)x(^+i))(tp). On the other hand, we allow v{j) to be equal to one or more of the 
v{t) when £ 7^ j ± 1. 



As outlined in the introduction, while in quantum statistical mechanics it is commonly 
expected that an expression as in (2.1) can always be reordered by bringing together 
operators at equal times, in this paper we would like to consider dynamical situations 
where the commutation relations between operators largely separated in time are of almost 
no use. Then, the natural algebraic structure to consider is that of a countable free- product 
2I00 = *i6No2li of copies of 21, which is the universal C*-algebra generated by an identity 
element I and by "words" w = X^\l\X^'^l,^ ■ ■ ■ X^"'\ that consist of concatenations of 

"letters" X^^'^ e 21. The subscript u^j) in ^^'(j) refers to which copy of 21 the letter X^^^ 
belongs. Concatenation, together with simplification rules, defines the product of words. 
More specifically, the rules for handling words are: for X, Y e % X E C, j e Nq and w, w' 
two generic words 

wljw' — ww' (2. 2. a) 

w{Xj + XYj)w' = wXjw' + XwYjw' (2.2.b) 
wXjYjw' = w{XY)jw' . (2.2.C) 



Notice that the product XY in (2.2.c) is not concatenation, but rather the usual operator 
product in the algebra 2t. Moreover, the adjoint w* of a word w = X^^^^X^^^^ ■ ■ ■ 
equals • • • 
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We shall refer to SIqo as to the asymptotic free algebra and equip it with an asymptotic 
state as follows. With the notation 

Avg{t ^ m) = ^lim - J2 m > (2-3) 

V any map from {1, 2, . . . , n} into Nq as above, we assume that the multiple average 
Avg(t ^ . . := Avg(t„ ^ • • • Avg(t2 ^ 

withn = 1, 2, exists. Then, we define a linear functional 0oo 

on 2loo by linearly extending 

the map on elementary words X^||^X^^2) ' ' ' -^^i4n)i 
</'oo(^i;i)4|)---4(y := Avg(t..</>(x«(^i))X(^)(t.(,))...X(-)(t.(„)))) . (2.4) 

The linear functional 0oo is such that it does not depend on the insertion of identities, 
namely 

Avg(t3 ^ Avg(t2 ^ Avg(ti ^ 0(1(^1)1(^2)^(^3)1^(^1)^(^3))))) = 
Avg(t3 ^ Avg(ti ^ cl>{x{U)Y{t^)Z{t^)yj) . 

Thus, (-^^3^1-^3) = 0OO (-^2 ^1^2) by relabelling the dummy summation indices. More- 
over, the linear functional </)oo is positive and one has the following general result: 

Proposition 1. The functional 4)oo defined in (2.4) extends to a state on 2loo- Further- 
more, 

0OO o ©s = 0OO and (poo <=> cte = (poo , 

where ©s satisfies &s{Xj) := (Qsj{X)) . for j G No and X G 21, whereas 6 is any order 
preserving infective transformation o/Nq andag is the *-homomorphism o/2too determined 
by ae{Xj) := ^^(j). 

We shall say that the asymptotic state 0oo is permutation invariant when 

ci>^(x%x^%---x^';\)=(l>^(x^^\,,x^^\,,---x^''\ ,) vx(^) e2l, (2.5) 

TOO I j^^i^ j^^2) / ^ \ 7roz/(l) 7rof(2) 7roi/(n) / ^ ^ ' \ / 

where tt : Nq i-^ Nq is any bijection. 
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3 Clustering properties and fluctuations 

The clustering properties of the dynamical system (2t, G, 4>) determine much of the struc- 
ture of the asymptotic state 0oo and hence of the statistics on the asymptotic free algebra 
2loo associated with (21, 0, <;^). 

Typically, two degrees of mixing are to be distinguished in quantum systems, called 
weak and strong clustering. They correspond to 

lim 0(xy(t)z) = 0(xz)0(y), X, y, zg21, (3.1) 

respectively 

(t){XY{t)Z S{t)T) = ^{X ZT)^{Y S), X,Y,Z,S,Te%. (3.2) 

t — >CXD 

The latter property, if it holds, is equivalent to hyper-clustering, that is to [NTl, ABDF] 

■ ■ -^'"'('"w)) = n<n.,„-.,,,^"") . (3.3) 

3 

where the limit is taken in such a way that all times and the differences between the 
different ones go to infinity. It is hyper-clustering that allows the reordering of multi- 
time correlation functions (2.1) with repeated times, when different times become largely 
separated. On the level of the asymptotic state (/)oo, strong (equivalently hyper-) clustering 
leads to: 

Proposition 2. Let (21, 0, (j)) he strongly clustering. Then, the asymptotic state 0oo defined 
on 2loo by (2.4) is permutation invariant. 

Essentially, if (3.2) holds, then using hyper-clustering, the order in which the single 
time-averages in (2.4) are performed does not matter for the averages coincide with the 
time-limits (3.3) of the multi-time correlation functions. Moreover, from o = it 
follows that (j)ao{Xj) = (j){X) for any j e Nq and X e 21, whence 

^oc(x:;> ) =n^(n.,„-,„^<") ■ (^■*) 

3 

The structure of the expectations on 't^^o calculated with respect to the asymptotic state 
0OO is compatible with imposing for j ^ k the commutation relations [2lj , 21a;] = on 2loo 
and corresponds to the usual notion of commutative independence of random variables. 

It is well-known that the notion of independence embodied in (3.4) is incompatible 
with that of free independence [VDN] which asks that the correlation functions of a state 
ip on a free product -kj^j of C*-algebras satisfy 

,/,(x<;'x£'...X<:')=0 (3.5) 
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whenever jk ^ jk+i and xj^^ is centred (i.e. ■0^-^)^'*^ = 0) for all k. 

The two notions of independence of above are somehow extreme. Many other possi- 
bilities for the structure of (poo may arise from the dynamical properties of (21, 0, (f)) and 
we shall try to expose some of them by looking at limits of the form 

hm </>^(F^(xW)F^r(X(2))...F^r(xW)) , (3.6) 
and establishing a central limit theorem for the local fluctuations Fn{X). 

Definition 1. Let N be a natural number and X e 21. v4 local fluctuation Fn{X) is the 
following element o/2loo 

Fr^iX) := ^ - 0(X)l) . (3.7) 

i=l 



By adapting an argument in [SvW], it can be shown that the cluster condition 

hm 0(z«(t,(i))..-ZW(t,(,))FZ(^+i)(t,(,.+i))---Z(-)(t,(„))) =0 , (3. 



inf \t^ — ^J 



for all centred y e 21 and v : {1,2, ...,n} i— Nq, which is stronger than weak cluster- 
ing (3.1), but weaker than strong clustering (3.2), is sufficient to ensure that only moments 
of even order contribute to the limit joint distribution of fluctuations. 

Proposition 3. Let us assume that (3.8) holds in (21,0,0). Then, with 0oo defined 
by (2.4) and X^^\ . . . , X^^^ in 21 centred observables, 

{0 r = 2n + 1 

;^I:''V«(X%...<£,) r = 2n. 

Y^^u^ means that we have to sum over all partitions v of {1,2, ... , 2n} into pairs ((cti. Pi), 
(0^2,^2), • • ■ , («n, /9n)) i-e. we choose sites aj < (5j such that i'{oLj) = I'iPj) = j with j 
running from 1 to n. 

As an immediate consequence, we also have 

Corollary 1. // the asymptotic state (poo is permutation invariant in the sense of (2.5), 
and X^'^\ . . . , X^'^^ are centred observables of%, then 

fO r = 2n + 1 



V, ord 



5 



where ^^^^ ' means that the sum is over all ordered pair partitions v = ((q!i,/?i), 
(ct2,/?2), • • • , (ctn, /^n)) o/ {1, 2, . . . , 2n}, 1.6. we choose sites aj < (3j such that I'icy.j) = 
I'iPj) = j with j running from 1 to n and a\ < a2 < ■ ■ • < OLn- 

Corollary 2. If the dynamical system (21, 0, 0) is strongly clustering, then the fluctuations 
Fn{X) of observables X e 21 such that (piX) = and 0(^^) = o^^^ tend to Gaussian random 
variables with zero mean and variance a. 

A crossing occurs in a given a pair partition v of {1, 2, . . . , 2n} into n pairs {aj, Pj), 
with aj < Pj, when 

aj < ak < Pj < Pk for some j, k E {1, 2, . . . , n} . (3.10) 

Denoting by c{i') the number of crossings in i/, we say that v is non-crossing if c{i') = 0. 
If is a non-crossing pair partition of {1,2,..., 2n}, then its pairs (ttj, Pj), j = 1, 2, . . . , n, 
are nested, that is if aj < < Pj for some j,k, then also Pk < Pj- According to 
expectations associated with pair partitions, one has various types of generalized Brownian 
motions [BKS,BS,S,SvW,vLM], that can be characterized by the contribution —l<q<l 
of each crossing 

n 

^441, ■ ■■Kpn)) = 1'"^' n ^^^'^'=^^^^'=0 • (3-11) 

k=i 

For q = ±1 the Brownian motion is termed Bosonic, respectively Fermionic, while for 
g = the g-deformed Brownian motion is called Free [VDN] and to it only non-crossing pair 
partitions do contribute. 

It turns out that if, besides the cluster condition (3.8), it also holds that, for any time 
independent choice of observables A and C and centred observables X, Y, B, 

Avg(t ^ (l)(^AX{t)BY{t)C^^ = , (3.12) 

then a free statistics for the fluctuations of temporal averages emerges. Namely, 

Proposition 4. Let (21, O, (/>) satisfy conditions (3.8) and (3.12). Then, with (poo the 
asymptotic state on 2loo defined by (2.4) and X^^\ . . . centred observables in 21, 

^hm </>^(f^(xW)F^(x(^))...F^(X(^-))) = ^E^'^ W) n H^^'^'^^^'^O ' 

where the sum extends over all non-crossing pair partitions v = (^{ai, Pi), {a2, P2), ■ • 
{an, Pn)) of {1, 2, . . . , 2n} and c{u) denotes the number of crossings in v. If the asymptotic 
state also permutation invariant, then 

^lim^</>oo(F..(xW)F^(X(^))...F^(X(^'^))) ^5o,c(.) H ' 

V k=l 
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where the sum now extends over all ordered and non-crossing pair partitions v o/{l, 2, . . 
2n}. 

Corollary 3. // the dynamical system (21, 0, (/>) satisfies conditions (3.8), (3.1) and the 
asymptotic state (poo is permutation invariant as in (2.5), then the fluctuations Fn{X) of 
observables X e 21 such that (f){X) = and = cr^ tend to semicircularly distributed 

random variables with zero mean and variance a. 



4 Quantized automorphisms of the torus 

We shall consider discrete time dynamical systems (21, 0, cj)) constructed as follows. Let u 
and V be two unitary operators satisfying the commutation relations 

uv = e^'""^ vu , 

where < ^ < 1 will be referred to as the deformation parameter. Then, we construct the 
Weyl unitaries 

indexed by two-dimensional vectors with integral components m = (mi, 777,2) £ Z^. The 
operators W^(m) satisfy the Weyl relations 

Wg{in)Wg{n) = e»'^^'^(™'") Weini + n) = e^*'^^'^^™'") W0{n)We{m) , (4.1) 

with a the symplectic form 

a{m, n) := min2 — ni2ni . 

Notice that VFe(m)* = Wo{—m) and that products njW^6i(mj) are reducible to just a 
single WffC^^ m^) multiplied by a phase. 

There exists a unique C*-norm on the linear span of the Weyl unitaries and the C*- 
algebra 21 is the completion with respect to that norm of {W0{f) \ f complex function on 
with bounded support}, where 

Weif):^ 5] f{ni)We{m) . (4.2) 



The dynamics O is determined by the linear extension to 21 of the map 

QiWeim)) -WeiTm) . (4.3) 

T = ( ^)isa2x2 matrix with integer entries such that ad — bc = 1 and a + d > 2. In 



c d 

this way, the eigenvalues A+ = A > 1 and A_ = A~-^ of T are irrational and the associated 
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eigenvectors v± correspond to stretching, respectively shrinking, directions with positive, 
respectively negative, Lyapounov exponent ± log A. Since cr(Tm, Tn) = (T(m, n), the map 
in (4.3) preserves the commutation relations (4.1) and thus extends to an automorphism 
of 21. 

Finally, the reference state is chosen to be the unique tracial state on 21 given by 

0(W^,(m)) :=So,m . (4.4) 

The state is clearly ©-invariant. 

Remark. Via the GNS construction an infinite-dimensional Hilbert space representation 
of 21 is obtained, even in the case of rational deformation parameters 6. This makes the 
quantization procedure introduced above drastically different from the ones presented in 
the literature [BBTV,D] where 6 is rational, typically 1/N and periodicity is imposed, 
namely 

In such cases, the resulting algebra acts on an A^-dimensional Hilbert space and the quan- 
tized hyperbolic automorphisms of the torus are a useful testing ground for studying the 
classical limit of quantized classically chaotic systems. Indeed, notice that by setting ^ = 
(or letting — > oo), the Weyl commutation relations (4.1) are implemented by the ex- 
ponential functions on the torus. This means that the algebra 21 becomes the Abelian 
C*-algebra of continuous functions on the torus which can be used to construct the corre- 
sponding algebraic classical dynamical system. 
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4.1 Number theoretical interlude 



Contrary to the finite- dimensional quantization procedures mentioned above, where quasi- 
periodicity spoils any true relaxation property, the dynamical systems (21, 0, (p) enjoy, 
depending on 9, sufficiently strong clustering properties. The behaviour of commutators 
of observables largely separated in time is determined by the value of the deformation 
parameter 



We (m) , We (T^n)] = (l - e^'^^^'^^'"'^'"^ ^We{m + T*n) , t 



e Z 



If 9 is rational, 9a {m, T*n) mod 1 can only assume a finite number of values, and hence no 
clear limiting behaviour occurs, unless of course ^ = 0. Some interesting behaviour occurs 
for peculiar irrational values of the deformation parameter. We shall first examine such a 
possibility in more detail and study the convergence mod 1 of the exponent 9a{m, T*n) 
for all m, n G Z^ when t ±oo. We adapt to the present problem a technique used 
in [AGL] and obtain as a sub-case a result derived in [N] with other means. 

Proposition 5. Let T = ( ^ ) with a,b,c,d^ Z such that ad — be = 1 and a + d > 2. 



c d 

For any m = (mi, 777,2) and n = (711, 712) £ Z^, define 

A(m, n) := (1 — a)(7ni77.2 + 77i7Tt2) — cmiTii — 6m2772 . 

Let A > 1 and he the irrational eigenvalues ofT, \± the corresponding eigenvectors 
and set 

/3(t) TrT* = A* + A"* , teZ. (4.5) 

Let cr(m, n) he the symplectic form m\n2 — m2ni on 7? x 7? and let 9 e [0, 1) he irrational, 
then, the limits 

Q±(m,n):= lim ^ f7(m, T*n) modi (4.6) 

t— »±oo ^ ' 

exist iff the following limit exists 

/?:= hm T2^/3(t) modi . (4.7) 
t— »-l-oo — 1 

When the limit in (4-.7) exists, its possible values are the rational numbers 

^'^Wr^' 0<r</3(l)-3. (4.8) 
In such cases, the limits (4-6) are of the form 

g!^(m, n) = ±l3r A(m, n) mod 1 . 
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Proof: Let v± be the normalized eigenvectors of T. They satisfy Tv± = A='=-'^v± with 
components with respect to the basis {(1,0), (0, 1)} given by 



v±{l) = h/^lP + {a-\±^Y , v±{2) = (A±i - a)/^h^ + (a - A±i)2 . 

By expanding m e as m = c+(m)v+ + c_(m)v_ and using (a — A) (a — A~-^) = —6c, 
we exphcitly compute 

/ \ min2A*"'"2 + m2niA* — A^+^fcmini + 6m2n2 + amin2 + am2ni) 
a(m, T n) = ^^^-^^ 

min2A~* + m2niA~*"''^ — A~*"''-^(cmini + 6777-2^2 + amin2 + am2ni) 

A^^l ■ 

Then, setting (f) -.^ 9 / {X^ - 1) , one works out that 

q+{in,n)— hm 0(min2/3(t + 2) + m2ni/3(t) — Ai(m, n)/3(t + 1) ) modi 

q-(in,n) = — hm 0( min2/3(— t) + m2ni/3(— t + 2) — Ai(m, n)/3(— t + 1) ) modi, 
t— >— 00 \ / 

where Ai(m, n) := a{min2 + 1211712) + crriini + 6m2?i2. Therefore, if the hmit in (4.7) 
exists, then the two limits above exist. Vice versa, if the previous two limits exist for all 
m, n e Z^, then, choosing proper values for m and n, the limit in (4.7) exists as well. 



The traces in (4.5) obey the recursion relations 

P{t)^P{l)p{t-l) - p{t-2), t>2, 

with P{0) = 2 and /?(1) = A + A~-^. Thus, if the limit in (4.7) exists, /? is determined by 

p - 2) = mod 1 , 

and ranges mod 1 among the rationals in (4.8). The explicit form of the limit q±{m, n) is 
obtained by inserting in the expressions for q±{in, n) the values of Pr- 



The previous proposition gives the explicit form of the limits in (4.7) if they exist. 
The next result concerns the values of 9 such that this is indeed the case. 

Proposition 6. Limt P{t)9 / {X^ — 1) mod 1 exists and equals (3^ = r/(/?(l) — 2), with 
< r < /3(1) - 3, iff 

^ = A£+(A-l)/3^ mod 1, £eZ. (4.9) 
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Proof: We define (j) := 9 / {\^ — 1) and divide the proof in three steps. 

First, we write (f> I3{t) = a{t) + 6(t), where a{t) is the largest natural number smaller 
than (p j3{t). If the Hmit in (4.7) exists, then Hmt 6(t) = for < r < /3(1) - 3. Let 
e{t) := h{t) — (3r and consider the quantities •0(t) := 4>P{t) — (3r — a{t) + e{t). They are 
easily showed to obey the recursion relations 

V'(t + 2) = /3(l)V'(t+l)-V(t) + r , t>0. 

The latter can be rewritten as 

a{t + 2) - f3(l)a{t + 1) + a{t) -r = f3{l)e{t + 1) - e{t + 2) - e{t) , 

whence we deduce that there must exist an integer T such that for t >T both sides of the 
above equality vanish. Indeed, the l.h.s. is an integer, while the r.h.s. goes to zero with 
t — > +00. In particular, this argument yields 

a{T + t + 2) = P{1) a{T + t+l)-a{T + t) + r , t>0 . (4.10) 

We prove by induction that for t > 0, 

t 

a{T + t + 2) := ^{t + 1) a(T + 1) - -f{t) a{T) + r ^ -i{k) , (4.11) 

fe=0 

with 7(t) the integers such that 7(0) = 1, 7(1) = (3{1) and 

7(t + 2) = /3(1) 7(t + 1) - 7(t) , t > . (4.12) 

The case t = is obvious. Suppose that (4.11) holds for < t < s — 1 and rewrite (4.10) 
with t = s as 



,(T + 5 + 2) = a(T + 1) (/3(l)7(s) - - 1)) - a{T) [mii^ - 1) - ^{s - 2)] 

s-2 
fc=0 



s — z 

T (l + /?(l) + ^(/3(l)7(fc + l)-7(fc) 



Using (4.12), a(T + s + 2) turns out to be of the form (4.11) with t — s. 
The second step consists in observing that the coefiicients 



c{t) 



A2-1 



of the expansion of (2;^ — fi{l) z + 1) ^ around 2; = fulfil the same recursion relations 
as the 7(t) in (4.12), with the same initial conditions. Therefore, ^{t) = c{t) for alH > 
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and one has 



m = (A^ - 1) 7(t - 2) + A-* (4.13) 

t-2 

7(t - 1) = 7(t - 2) + 1 + m) - 2) ^ 7(A;) t>2 (4.14) 

fe=0 

A-* = 7(t) - A 7(t - 1) , t>l (4.15) 



We can use the previous relations to prove that the hmit in (4.7) exists and equals /3r 
in (4.8) iff 



A a{T + 1) - a(T) + /?r(A - 1) mod 1 . (4.17) 



We write 4> = / (A^ — 1) and (j)[3{T + t) = a{T + 1) + h{t) as in the beginning of this proof. 
Then, dividing by /9(T + t) and letting t +oo, one recovers (4.17) by means of (4.5), 
(4.11), (4.13) and (4.16). Vice versa, if 9 is as in (4.17), (4.5) and (4.13) yield 

A 



lim (3{t) (x a{T + 1) - a{T)] mod 1 = 0, 

t^ + CX3 A — 1 V / 



lim I3{t) A"^ I3JX - 1) mod 1 = lim I3ril3{t + 1) - (5(t)) mod 1 

t— »+oo t^+oo ^ ' 

= lim /3^(7(t- 1) -7(t-2)) modi, 

t— >+oo ^ ' 

so sufficiency follows from (4.14), since Pr = r — 2). 

In the third and last step we reduce the expression (4.17) to the simpler form (4.9) by 
using the fact that, in (4.17), the integers T, a{T + 1) and a(T) are so far undetermined. 
In fact, the sufficiency of (4.9) to guarantee the existence of the limit (4.7) is proved along 
the same lines as for that of (4.17) above. For the necessity, we have to show that (4.17) 
reduces to (4.9). Using (4.15), (4.17) reads 
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= A d + /3^(A - 1) (^7(T) - A 7(T - 1)) mod 1 , 



where d := 7(T - 1) a(T) - 7(T - 2) a(l + T). By using (4.14) and the fact that 



we can rewrite 



(A-1)(7(T)-A7(T-1)) =/3,(A-l) + di + A5, 
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where di := d + r Ylp^Q l{p) and 5 := —r Yl^=o l{p)- Since d\ and 5 are generic integers, 
the expression (4.9) is obtained with £ = d + di. 



The following and last number theoretical result is a simple characterization of certain 
powers of the evolution matrix T. 

Proposition 7. Let T be a 2 x 2 hyperbolic matrix with integer entries, unit determinant 
and trace > 2 as in Propositions 5 and 6. Then, the matrix T^(i)~2 is congruent to 
the identity matrix mod — 2), in the sense that its diagonal entries are congruent to 
1 mod — 2) and the off-diagonal ones are congruent to mod — 2). 

Proof: The eigenvalues A and of T solve — I3{l)z +1 = 0, therefore they satisfy 
the matricial recursion relations 

= /^(l) n x-(n-i) - n x-(n-2) for n>3. 



A- ; - \ A-(-i) ; ~ V ^"^""'^ , 

Because the eigenvalues of T are non-degenerate we can find a similarity transformation 
which allows us to diagonalize T, such that the above equality may be turned into 

rjrn /3(l) rjTn — 2 

Then, defining 5n := - I 

Sn = Pil) Sn-1 - Sn-2 + - 2 , with = and Si = T-l, 
the differences Sk — Sk-i, k > 2, are connected by the following relation 

Sk - Sk-i = Sk-i - Sk-2 + {m - 2) {Sk-i + 1) , 
whence, by telescopic summation, 

n 

5,= (/?(l)-2)(J]5fc_i + n-l)+5,_i + 5i . 

k=2 

Since all matrices Sn have integral entries, Sn = Si+Sn-i mod (/9(1) — 2) and, by iterating 
n times the previous congruence, 

Sn = nSi mod(/3(l) -2) . 

Thus, when n = /3(1) — 2, Sn turns out to be congruent to mod {/3{1) — 2) and the result 
follows. 
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4.2 Clustering properties and statistics of fluctuations 



We return now to the discussion of the ergodic properties of the quantized hyperbohc au- 
tomorphisms of the torus (21, Q, (j)) and of their dependence on the deformation parameter 

^e[o,i). 



Proposition 8. For all 6 e [0,1), the quantized hyperbolic automorphisms of the torus 
(21,0,0) are weakly clustering in the sense of (3.1). They are strongly clustering in the 
sense of (3.2) iff, according to Proposition 5, 

hm ^_^/3(t) = modi , 
that is, according to Proposition 6, iff = X£ mod 1, with £ e Z. 

Proof: The state (f) is tracial and the operators Wg{f) in (4.2), with finitely supported / 
on Z^, are uniformly dense in 21. Therefore, by linearity, the cluster property in (3.1) is 
proved by showing that 

lim (l>(w0{m)W0{T'n)W0{p)) = 

„7ri6lfCT(m,T*n)+CT(m+T*n,p)) r n. 

hm e V ^ ' ^ V do^+T^n+p = , 

t—^00 

for generic integer vectors m, n and p e Z^. The commutation relations (4.1) and the 
definition (4.4) of the state have been used to derive the equality of above. The limit 

is equal to zero because of the hyperbolic character of the matrix T. Indeed, by means of 
the eigenvectors v±, one puts into evidence the Lyapounov exponent log A, 

T*n = c+(n)A*v+ + c_(n)A-V_ . 

It is thus evident that in the limit of large \t\, the condition m + T*n + p = cannot be 
fulfilled. 



Analogously, one calculates 



Weim) , W0{T*n) * W0{ia) , W0{T*n) ) = 4sin2(7r^(7(m,T*n)) 



Strong clustering (3.2) implies that the right hand side of the above equation must vanish, 
that is that lim^ ^cr(m, T*n) mod 1 = 0, for all m and n G Z^. On the other hand, if 
lim^ 9a {m, T*n) mod 1 equals for all integer vectors m and n, then all commutators as 
the one above vanish in norm and therefore the dynamical system is even norm-asymptotic 
Abelian and hence, as it is weakly clustering, it is also strongly clustering. 
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Next we shall compute the asymptotic state (pco and the fluctuations for the hyperbolic 
toral automorphisms (21, 0, 0) with deformation parameter 9. We shall sec that, depending 
on the choice of the deformation parameter, these aspects can be quite difi'erent. We shall 
rely on formula (2.4) with single-time averages as in (2.3). The "words" in the asymptotic 
free algebra consist now of free products of elements of the type W6i(/*^-'^),y(j) with W0{f) 
given by (4.2). Because of linearity of / i-^ W^e(/) we may restrict ourselves to studying 
multi-time correlation functions of Weyl operators 

t ^ 0(|M^e(T*"'<''nW) Vre(T*-'(2)n(2)) . • • (T*-'(") n^^) 

where v maps {1, 2, . . . , n} into {1, 2, . . . , s}. 

The following lemma can be found in [NT2] . For the sake of completeness we provide 
the proof as well. 

Lemma 1. Let T be any hyperbolic matrix defining the dynamics of (21, ©,</)) as in (4-3) 

and let a be the symplectic form a{m, n) = min2 — m2?^i with m, n G Z^. For almost all 
9 G [0, 1), the orbits /c G Z 9a{m,T'^n) are uniformly distributed over the circle T of 
perimeter 1. 



Proof: Uniform distribution means that V/ G C(T) we get 

i ^ 
m — 

N- 



1 ^ 

k=l 



This is equivalent to requiring that 

1 ^ 

fN{9):^^Yl e''^^^^'^^"^'^'") ^0 Vn G Zo . 
fc=i 

To see whether this holds for almost all 9, we consider 

/ d9\JM' = ^, E / d9e'-^-''^("^'^^'-^'^-K 
Because cT(m, (T^ - T^')n) 

is an integer, we have that 

^Q^2nin9a{m,{T''-T'^')n) ^ f 1 if f7(m, {T'' - T''')n) = 

Jo 1 else 

And, because a{ni, [T'' - T''')n) ^ y m e Z'^ ^ T'^n = T''' n ^ k = k' , we have that 

AT 

1/N — > for iV ^ oo . 



fd9\fM\' = W E 1 = 
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Using the previous lemma and Propositions 5 and 6, we deduce 



Corollary 4. a) There exists a set Z C [0, 1) o/ measure 1 such that for all 9 ^ Z and 
for all m, n ^7? Q), we have 



1 ^ 



h) //lim^ ^cr(m, T*n) mod 1 exists, then 

T 



lim LY^'^-^Oo{r^,T^r.) 



1 

lim — "V" g27riea(m,T*n) ^ 27ri/3^A(m,n) 



where f3r and A(m, n) are defined in Proposition 5. 

In order to proceed with the construction of an asymptotic state on 2too following the 
prescription of (2.4), we need a preliminary technical result which is proved in [ABDF]. 

Lemma 2. For d,k E J>i define 

Ak(+ + \ _ / if \ti — tj\ < d for some 1 < i ^ j < k 

Then, if the multiple average (2.14) of a uniformly bounded function f : — > C exists, 
we have 

Ayg(tk ^ •••Avg(t2 ^ Ayg(ti ^ /(ti, ^2, • • • , tfc))) ■ ■ ■) = 

Ayg(tk ^ •••Avg(^t2 ^ Avg(^ti ^ A^(ti,...,tfc)/(ti,t2,...,tfe))) •••) • 

Proposition 9 and the subsequent corollary deal with the asymptotic structure emerg- 
ing from the first possibility of corollary 4, while proposition 10 and its corollary deal with 
the second case. 

Proposition 9. Let us consider a hyperbolic toral automorphism determined by a defor- 
mation parameter 6 E Z as in corollary 4, po^rt a). The multiple average in (2.4), with 
single-time averages as in (2.3), exists and the state 0oo it defines on 2loo permutation 
invariant as in (2.5). 

Proof: Because of linearity, we need only to consider expectations of words of the type 
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which are to be computed as the multiple averages 



AYg(ts ^ ■ • • Avg^t2 ^ Avg^ti 1-^ 

V maps {1, 2, . . . , n} into {1, 2, . . . , s} and we may assume that all Weyl operators are 
different from the identity, that is that none of the n^-?) = 0. 

Before taking the average we group all Weyl operators belonging to a same time, gathering 
hereby a phase factor according to the commutation relations (4.1). Then, using Lemma 2, 
we reduce the problem to the computation of the multiple time-average of the phase factors. 
By the assumption on the deformation parameter ^, these averages are always zero unless 
they are time- independent. This proves the permutation invariance of 0oo- 

With p G {1, 2, . . . , s}, let Ip denote the set of natural numbers i G {1, 2, . . . , n} such that 
— p. For any given p G {1, 2, . . . , s}, the Weyl operators We(T'*'''''n*^*)), with i G /p, 
can be brought together, from right to left, to form the composite word 

Using (4.1), the regrouping produces a phase factor 

Help k€K{i) jeJi{k) ^ 

The index set K{i) specifies those times 7^ tp that are encountered while commuting 

over the various W6i(T*'^(*)n*^*-') in order to concatenate the former 
word to the right of a previous word VF6((r*pn*^-^)), j G Ip. 

We shall regroup the Weyl operators by first bringing together all those at time ti, then 
all those at time ti and so on. We end up with 

F(ti,t2,...,t.)</>(w^/i(ti)H^/,(t2) ■■■WiXts)) , (4.18) 

where 

s 

F(ti,t2,...,ts^ = Yi Fp{tp,tp+i,...,ts) with 
F,(«p,tp+i,...,t.) = exp( J] Yl Yl 2iie<7(T"nO),TW")j , 



Help keKp+i{i) jeJi(k) 
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and Kp+i{i) labels times tk with k >p+l. Next, using Lemma 2, we can always restrict 
ourselves in the multiple timc-avcrages to times t,y(j) that are sufficiently separated from 
one another. As the dynamics n i-^ Tn is hyperbolic on Z^, the separation between 
different times can always be chosen in such a way that 

u{i) ^ =^ T*'^(i)n'^i) + T^-^Wn'^^) ^ o and 

n 

^T*-«)n(-?) = ^ Hp := J]nW = = 1, 2, . . . , s . (4.20) 

j=i ieip 

Prom (4.20) and (4.4), we see that the only possibly non- vanishing averages are those for 
which Hp = for each p separately, in which case Wj^ = Gp Wg{np) = Gp I, Gp being a 
suitable phase obtained through (4.1). We can therefore write 

s 

(n G<i^o,n,) Avg(t, ^ • • • Avg(t2 ^ Avg(ti ^ t2, . . . , t.))) ■ ■ ■) . (4.21) 

We now take the average with respect to ti. In the product (3.19), only the factor Fi can 
depend on ti . Since 9 belongs to Z, this average is either zero or one. Moreover, it can only 
be one if Fi does not depend on ti . We can now successively average over the consecutive 
times t2,t3, ■ ■ ■ and conclude that either the average is zero or that it is independent of all 
ti,t2, ■ ■ ■ ,ts. In both cases, 0oo is permutation invariant in the sense of (2.5). 



Corollary 5. A hyperbolic toral autom,orphism determined by a deform,ation parameter 
9 & Z as in Corollary 4, part a), satisfies condition (3.8) and (3.12). Therefore, the fluc- 
tuations of a centred self-adjoint element WQ{f) in (4-2), are semicircularly distributed. 

Proof: Condition (3.8) is satisfied because of the hyperbolic character of the matrix T 
which implements the dynamics. Thus, because of linearity and since the reference state 
is the trace, it suffices to check condition (3.12) in the form 

Avg(t ^ <j){We{p)We{T'ni)We{n)We{-T'in)'^'^ = , (4.22) 

for all m, n, p e Z^. In fact, 

a) products of Weyl operators can be reduced to a single Weyl operator multiplied by a 

phase factor; 

b) Weyl operators on the right of Wg{—T*m.) can be moved to the left of Wg{T^m) by the 
tracial properties of the state (f) and 

c) the two and only two Weyl operators at time t carry opposite integral vectors because, 
otherwise, the hyperbolic character of the matrix T and (4.4) would force the expectation 
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and thus also its average, to vanish asymptoticaUy. 
Then (4.22) is obviously true, as 

Avg(t ^ (j)(w0{p)We{T'm)W0{-n)Wei-T'm)'^^ = 

(f)(w0(p)W0{-n)^ Avg(^t^exp(^27rz^(7(T*m,n))) . 

On the basis of this result, we want to compute limjv 0oo ^-Piv(W6i(/))^'^^ where Fat is the 

local fluctuation defined in (3.7). Notice that W0{f) is self-adjoint iff, for any n e Z^, 
/(n) = /(— n) and it is centred iff /(O) = 0. We shall denote by the support of /. Since 
is permutation invariant, we use Corollary 3 and obtain a semicircular distribution for 
W0{f) with variance 

^ = <^(w^^(/f) = ll/lli = El/(q)l'- 



Remark The previous proof could suggest that the state 0oo which was constructed in 
Proposition 9 is a free product of traces, namely that 

^.(4"xf...x<:')=o 

whenever ^^xj^-'j = and jk ^ jk+i for all k. This is, however, not the case. An easy 
counterexample is obtained by considering a correlation function that is independent of a 
time appearing in the product of observables such as 

t ^ (j)(w0{n)W0{T^n)W0{-2n)W0{T*n)W0{n)W0{-2T*n)^ n ^ . 

Each of the observables in the correlation function is centred but, instead of vanishing, 

(f>oo[W0{n)iW0{n)2W0{-2n)iW0{n)2W0{n)iW0{-2n)2) = 1 . 



We shall now see that, when the deformation parameter 9 is chosen not to belong to 
the dense set Z, but, instead, equals any of the special values for which the second part 
of Corollary 4 holds, then the multiple time- average (2.4) still defines an asymptotic state 

but, except when = 0, it is not permutation invariant. 

Proposition 10. // the conditions of Corollary 4, P^^^i b), are fulfilled, then the multiple 
average (2.4) exists and it defines an asymptotic state 0oo on the asymptotic free algebra 
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Proof: We can follow the proof of Proposition 9. The asymptotic state 0oo is determined 
by the expectations 

,A„<.(w.(„<'>)„„,w',(„<^>)„,„...w.(„«">)^,„,), 

where v maps {1, 2, . . . , n} into {1, 2, . . . , s}, s < n, and the vectors n*^-^^ G T? may be 
supposed ^ 0. By suitable regrouping as in formula (4.18), the expectations above are 
well-defined if the multiple time-average (4.21) exists. Thus, we start considering the time 
t\. By virtue of the assumption on the existence of the limits of the exponents, we get 

Avg(ti ^F(ti,t2,...,t.)) = 

AvgTti ^exp^'j^ Yl 27rz^f^(T*'=n(^), 

P=2 ^ Hell keK2ii) j€Ji(k) 

s 



llFp{tp,...,ts) expfYl E Y ^'''f^rA(T'^n^^\n^^^)Y (4.23) 
p=2 \eiikeK'2{i)jeJi{k) ^ 

In the last equality, we have used the quantities (ir and A(m, n) which have been introduced 
in Propositions 5 and 6. Furthermore, the index sets Ki{i\ i E Iq and £ > q + I contain 

k>e. 

When averaging with respect to t2, the factor F2{t2,ts, . . . ,ts) tends asymptotically to 
a term similar to the exponential in formula (4.23) of above, without t2 dependence. 
Therefore, we only have to compute the average with respect to t2 of a contribution of the 
form 



exp 

T. Ar- 7. /'0\ / 



which may come from the presence of the time t2 among those indexed by A; e ^2(1) in 
the exponent in (4.23). 

We rewrite the quantity A(m, n) introduced in Proposition 5 as the scalar product A(m, n) = 
(m, Sn) where 5 = ^ ^ ^ b"^^ Since powers of T transform into and S has 

integral entries, we use Proposition 7 to deduce that 

^(j.m{m-'2)+s^U) ^ n«) = {T'n^^\ 5n«) + (/3(1) - 2) iV(n«, n^^')) , 

for all m e Z, with < s < /3(1) - 3 and iV(nW,n(^)) a suitable integer. Moreover, 
(3r — r I (/3(1) — 2) with < r < /3(1) — 3, thus the following cyclic properties hold 

^^p(E E 27rz/?,A(T-(^(i)-2)+-n(^),n«)') = 

^^'7^^^^^^ , (4.24) 

expf^ E 27rz/9^A(^r«n(^), n(^)) 
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After the first two averages, we remain with 

P=3 ^i€l2k€K3{i)jeJi{k) ' 

^^p(E E E 27ri/?,A(T*^n(^),nW)) I)? 



Se/i keK3{i) j£Ji{k) 
where, using (4.24) 



Dl := Avg(t2 ^ expT^ ^ 27rz/?^A(^T*^n(^), n^^^)") 

If we go on and consider the average with respect to ts, the result is 
Avg^tg ^ Avg(t2 ^ Avg^ti ^F(^ti,t2,...,ts)))) = 

nFp(tp,...,t,) expl'j] E 27ri/?,A(T*'=n(^),n«)] 

^""pfE E E 27rz/?^A(|T*'=n(j), n^)") x 

\el2k€K4(i)j€Ji{k) ^ 

^^p(E E E 27rz/?.A(T*'=n(^),n«))i?Li;?, 



X 



where 



1 '^^^^"^ / r 

^'^■■=3(^2 E exp(2.^/5J^ A(T^nO),n(^)) + 

J2 A(TV'^),n(^))y 



After averaging with respect to the remaining times ts, ^4, . . . , ts, we obtain a well-defined 
asymptotic state 0oo given by 

.A„o(»'.(n<'l)„,„iy4nP')^p,... = n 

p=2 

where 
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The asymptotic states just constructed are not permutation invariant as in (2.5). 
Indeed, when ^i 



W0{p)3 We{p)i W0{m)3 W0{-m)2 We{-p)i W0{ra)2 Wei-p - 111)3 j 7^ 
0OO [Weip)^ We{p)3 Weini)^ W^e(-m)2 We{-p)3 Weimh Wei-m - p)i) . 

By relabelhng the summation indices, as in the proof of Proposition 9, the inequahty of 
above comes about if we show that, with shortened notation for the successive averages. 



AvgsAvg2Avg^i^<P(We{T*^p) We{T*^p) We{T'^ni) W0{-T*^in) W0{-T*^p) 

We{T'^m)We{-T'^{m + p))'^ ^ 
Avgi Avg2 Avg3((/>(vre(T*^p) We{T''p) W^^T^^m) 

Wei-T'^m) Wei-T'^p) We{T'^in) We{-T''{in + p))^ . 

We now bring together from right to left the words belonging to ti, t2 and t^, thus obtaining 
We (T*3 p) Wo (T*^ p) We (T*^ m) We (-T*^ m) We {-T'^ p) We (T*^ m) We (-T*^ (m + p)) = 
exp 

Then, arguing as in the proof of Proposition 10, we compute 

Avg3Avg2Avgi(^exp(^27rz^(a(T*^m,T*ip) -a(T*3in,T*ip))^^ = 
Avg3Avg2(^exp(^27ri/3^(A(T*^m,p) - A(T*3m,p))^^ = 

/3(l)-3 

cxp(27vi(3r^{T'm,p) 



1 



/?(!)- 2 



On the other hand, 



s=0 



(4.25) 



Avgi Avg2 Avg3 (^exp (2^6 (a (T*^ m, T*i p) - a (T*^ m, T*i p) ) ^ ^ = 



Avgi Avg2 l^exp (2m0a (T*^ m, T*i p) j exp \^2T:iPr ( A (T*i p, m) ) 
Avgi ( exp (^-27rz/5r A (T*i p, m) ) exp (27ri/3rA {T*^ P, m) H =1 



Clearly, the expression (4.25) can be made different from 1 by suitably choosing m and p. 
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Notice that, if (3r = 0, which is the case when r = 0, then (4.25) equals 1 as well. 
In fact, Proposition 8 tells us that the system (21, G, (p) is strongly asymptotically Abelian 
and Proposition 2 ensures that the asymptotic state is then automatically permutation 
invariant. 

In order to characterize the statistics of fluctuations in the asymptotic algebra 2loo) 
when the asymptotic state (/)oo is defined as in Proposition 10, we show that, besides being 
weakly clustering, the dynamical systems (21, Q, (j)) satisfy condition (3.8). 

Proposition 11. Let (21,0,^) be a quantized toral automorphism with deformation pa- 
rameter $ such that Unit 9(3 {t)/{X'^ — 1) mod 1 = f3r as in Proposition 5, with r ^ 0. 
Then, 

. lim </>(w,(/W)(t.(i)) • ■■Woif^^^)it,a))Wo{9) 

mt \ti—tj |— >cxD \ 

W^e(/^^ + '^)(tKi+l))---^^(/^"^)(^^(n))) = , 

where W0{f^^^), W0{g) are defined as in (4-2), by means of finitely supported functions, 
and 

We{f^'^){K^^)) = E /(^)(m)I^,(T*^(^)m)) . 

Proof: Because of linearity and the finiteness of the supports of /^^\ we can restrict 
our considerations to expectations of the form 

V^A^(^T*-'(.+i)p(ij+i)j ...Vt7^(^T*-(")p("))^ =0 , 

where m 7^ 0, for the corresponding word is assumed to be centred. If we now group 
together Weyl operators at equal times, T*m is not matched by any of the other operators. 
On the other hand, using Lemma 2 and considering the smallest difference — 
between different times sufficiently large, we can always force upon the supports of Weyl 
operators at equal times a condition as in (4.20). Otherwise, the expectation would vanish. 
But then m = 0, which is impossible. 



The result of above shows that condition (3.8) may hold in systems which arc weakly, 
but not strongly clustering. Therefore, quantized toral automorphisms (21, 0, (f)) with de- 
formation parameter 9 ^ Z (see Corollary 4) fall in the class of dynamical systems whose 
time-asymptotic fluctuations can be handled via Proposition 3, by means of pair partitions 
only. Condition (3.12) obviously does not hold for any (3r since its equivalent version (4.22) 
is easily violated. Taking into account that, when (3r 7^ 0, the asymptotic states 0oo are not 
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permutation invariant, we can summarize our findings for the previous class of quantized 
automorphisms of the torus in 

Corollary 6. When the deformation parameter 9 is such that lim^ /3(t)6'/(A^ — 1) mod 1 
equals (5r with (5r as in Proposition 5, the statistics of time- asymptotic fluctuations of 
quantized hyperbolic automorphisms of the torus (21,0,0) is as follows. For (3r = 0, 
the fluctuations (3.7) of centred observables We{f) are Gaussian random variables with 
variance a = H/Hl- For /3r 7^ 0, they obey a distribution law with vanishing odd moments 
and even moments M2n given by 

M2n :=lim0oo(i^iv(We(/)'")) = ^$]^'Voo(w^e(/)Kl)W^e(/).(2)---W^e(/).(2n)) 

where the sum is over all pair partitions of {1,2, ... ,n}. 

Proof: The statement follows from Propositions 3 and 11 and from Corollary 5. 
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